Universal monopole scaling near transitions from the Coulomb phase 
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Certain frustrated systems, including spin ice and dimer models, exhibit a Coulomb phase at low temperatures, 
with power-law correlations and fractionalized monopole excitations. Transitions out of this phase, at which the 
effective gauge theory becomes confining, provide examples of unconventional criticality. This work studies 
the behavior at nonzero monopole density near such transitions, using scaling theory to arrive at universal 
expressions for the crossover phenomena. For a particular transition in spin ice, quantitative predictions are 
made through a duality mapping to the XY model, and confirmed using Monte Carlo simulations. 
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In frustrated systems, where ordering is hindered by com- 
peting interactions, large fluctuations can persist even at tem- 
peratures T where the degrees of freedom are strongly corre- 
lated. This can lead to the formation of so-called "spin liq- 
uids" yj], characterized by topological order and fractional- 
ized excitations. These exotic features also have important 
implications for transitions out of the spin liquid into conven- 
tional phases. 

In particular, continuous transitions from a spin liquid have 
critical properties that are not given by a standard Landau de- 
scription (|2[] in terms of long-wavelength modes of an order 
parameter Currently known examples include columnar or- 
dering in the cubic dimer model t3,-i8jL the Kasteleyn transi- 
tion of spin ice ^^ in a (100) field lflll[l2ll, and Neel order- 
ing of the Heisenberg model on the pyrochlore lattice 11311. A 
recent study of spin ice 11411 proposed a set of such transitions 
induced by application of appropriate perturbations. 

In these examples, the spin liquid occurring above the crit- 
ical temperature Tc is described by the Coulomb phase of an 
effective U(l) gauge theory jlSll . This phase, with character- 
istic dipolar correlations, occurs when the system is restricted 
to a low-energy manifold satisfying a constraint on the diver- 
gence of an appropriately defined degree of freedom. The 
transition results from a perturbation acting within this mani- 
fold, of magnitude V, with Tc « V. There may be order for 
T < Tc, in the sense of spontaneously broken physical sym- 
metries (e.g., columnar ordering in the cubic dimer model), 
but this is not required, either for the existence of a phase tran- 
sition or for the present analysis. 

The two phases are distinguished by their response to the 
introduction of defects in the divergence constraint, inter- 
preted as charges, or "monopoles", in the effective gauge the- 
ory. A pair of opposite sign introduced into the defect-free 
system feels a long-range interaction of entropic origin. In 
the Coulomb phase (in three dimensions, 3D), this takes the 
characteristic 1/r form for large separation r, while it grows 
(at least) linearly with r, signaling confinement, in the low- 
temperature phase. (In the case of spin ice, the gauge charges 
are also physical magnetic monopoles lll6ll .) 

While a test pair allows one to distinguish the phases, any 
finite density of monopoles screens the dipolar correlations of 
the Coulomb phase. In other words, if the energy cost for a 
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FIG. 1: Schematic phase diagram for a system with a continuous 
transition out of the Coulomb phase. The latter is stable only at z = 
(blue line), so the confinement transition is an isolated point in the 
phase diagram. It nonetheless influences properties in a broad region 
and leads to universal scaling forms such as Eq. (|2). The dotted 
line (T > Tc) indicates a crossover from Coulomb-like behavior to 
a conventional paramagnet, while the dash-dotted line (T < Tc) is 
either a transition (in a conventional universality class) or a crossover. 
Both have the form z ~ \T - Tcf, where is a crossover exponent. 
The (red) aiTow is an example path as T is reduced at fixed A and V. 



monopole is A, the U(l) spin liquid exists as a distinct phase 
(in 3D) only for fugacity z = e"^^^ = JITII . Thermal phase 
transitions from the Coulomb phase therefore appear to in- 
volve a contradiction: nonzero temperature is required for the 
transition, but for any T > thermally excited monopoles 
render the Coulomb phase unstable lilSil . 

This work shows that, on the contrary, nonzero monopole 
density provides a valuable new perspective on the unconven- 
tional critical behavior. This can be appreciated by consider- 
ing z and T as independent parameters, as in Fig.[T] While the 
critical point exists precisely at T = Tc and z = 0, it has con- 
sequences for behavior in a region surrounding this point. In 
particular, scaling theory strongly constrains the dependence 
of physical observables on z and T, in ways determined by 
the properties of the zero-monopole critical point itself. The 
analysis at z > is close in spirit to the study of critical phe- 
nomena in systems of finite size L, where L"' = at the crit- 
ical point lie 



and of quantum criticality at T > II2II1 . Our 
universal predictions can be tested in numerical simulations 



such as those of the cubic dimer model ||3|], while understand- 
ing the behavior at z > is crucial for realizing this class of 
unconventional critic ality in experiment. 

In addition, this work provides numerical evidence for the 
existence of a new continuous transition in this class, occur- 
ring in a model of s pin ice and previously predicted using an- 
alytical arguments lll4ll . This transition has the advantage that 
its critical exponents are known to high accuracy, due to a du- 
ality mapping to the XY model. It is therefore possible to 
provide quantitative predictions for the expected critical be- 
havior, including at nonzero monopole fugacity, and these are 
confirmed using Monte Carlo (MC) simulations. 

We consider thermal phase transitions (i.e., T > 0; the 
model may be classical or quantum mechanical) in 3D be- 
tween a spin liquid and a "conventional phase". Specifically, 
suppose that for T > Tq and z = 0, the long-wavelength de- 
scription is a U(l) gauge theory in its Coulomb phase, while 
for T < Tc t he g auge theory is confining. Some transitions of 



this type fld. ll4ll can be viewed as Higgs transitions 1221123 1 
driven by the condensation of an emergent matter field (dual 
to the monopoles), while critical theories for others are found 
through mappings to effective quantum models 141 IllL 11211 . 
The precise nature of the critical theory is not of consequence 
for this analysis, which assumes only that such a theory exists. 

The phase structure for z > 0, illustrated in Fig.lT] depends 
on the nature of the confined phase at T < Tc and z = 0. If 
this breaks no symmetries, then it connects smoothly to the 
paramagnet. Otherwise, a phase boundary separates the para- 
magnet at r > Tciz) from the ordered state at T < Tc{z)- The 
latter transition belongs in a "conventional" universality class 
(or is of first order), in contrast to the point at z = 0. 

In terms of renormalization-group (RG) theory, the insta- 
bility of the confinement transition implies that turning on a 
nonzero z amounts to adding a relevant perturbation. There 
therefore exists a "scaling field" z, that is an increasing func- 
tion of z vanishing at z = 0, and that is conjugate to an eigen- 
operator of the RG transformation [20]. In other words, near 
the fixed point at f = 0, z = 0, rescaling by a factor b replaces 
the effective values of t and z by tV' and zV~. As usual y, > 0; 
the relevance of monopoles implies that y^ > 0. 

At most conventional fixed points, a perturbation appears 
as an additional term in the Hamiltonian, and symmetry often 
dictates that the scaling field be proportional (to leading order) 
to its coefficient. Here, nonzero z instead reduces to a finite 
value the energy cost of a monopole, oc |In z|, and there is no 
associated symmetry. That the appropriate scaling field is in 
fact simply z is therefore not obvious, but will be estabUshed 
using a mapping to a conventional ordering transition ll24ll . 

For concreteness, consider a classical model with discrete 
degrees of freedom Bf defined on the links ^ of a lattice 11511 . 
and a partition function 
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where the action S depends on the full set [Be]. The integral 



over 9i e [-n,n) constrains «, = div,B, where div,- is the 
lattice divergence at site /. 

First, consider the case where z = 0, so n, vanishes. Us- 
ing the identity 2i^/div,B = - Y^t Be grad/. 6, where grad^ 
denotes the lattice gradient on link {, one can in principle 
perform the sum over B(. Crucially, while the effective ac- 
tion S^^{d) depends on S, it can always be written in terms 
of gradf 0, so has an XY symmetry under global shifts of 0,. 
For appropriate <S, the angle variables can therefore order, cor- 
responding to the Coulomb phase, as in the standard duality 
between the X Y model and the current-loop model ll25ll . (As 
usual, the high- and low-temperature phases are exchanged 
by this duality mapping.) For a general interaction S, this or- 
dering transition need not be in the same universality class as 
the XY model (and the ordered state need not have uniform 
(e'*'), for example), but it remains true that the XY symmetry 
is spontaneously broken at the transition. 

For z > 0, the sum over n,- gives an additional contribution 
to S'^^iff) that breaks the XY symmetry explicitly and elimi- 
nates the possibility of an ordered state. To leading order in z, 
the result is a term -2z 2; cos 0,, corresponding to an applied 
field lixY f^ z acting on the angle variables. This mapping 
has therefore related the unconventional critical properties to 
a standard example of crossover scaling ll20ll . completing the 
argument that z is the appropriate scaling field. 

Scaling theory then implies that, for example, the singular 
part of the reduced (i.e., divided by temperature) free-energy 
density /s obeys 
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where t - (T-Tc)/Tc is the reduced temperature, a - 2-d/y, 
is the specific-heat exponent, and (p = yz/yt is a "crossover ex- 
ponent". (The a priori unknown function Oj- also depends on 
the sign of t.) Eq. ^ implies that nonzero z has greatest effect 
when |f 1"^ < z, while otherwise the argument of 0+ is small and 
/s may be approximated by its z = behavior Equivalently, 
one can define the "monopole screening length" /},„ ~ z ''^'*, 
where v is the correlation-length exponent; nonzero z becomes 
significant when the correlation length ^ ~ \t\^^ exceeds /!,„. 

The exponent (p, though governing the behavior at z > 0, is 
a property of the fixed point itself, and can be determined ex- 
actly at z = using an appropriately chosen correlation func- 
tion. At a fixed point, the correlations of a scaling operator 
(fi obey {ip(r)ip(r')) ~ jr - r'l"-*"'"-''** 12011 for large separation, 
where d is the spatial dimension. The correlation function 
corresponding to the monopole fugacity z can be identified by 
replacing z in Eq. ([T]i with a nonuniform source z,, 
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taking derivatives with respect to z, and Zj, and setting all 
sources equal to z. The resulting correlation function. 
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becomes, at the fixed point, G,j(0) = [-3(0)]"' Yj'\b i ^ '^- ^^^ 
sum is over configurations with |div,' B\ = 5,',- + 6,' j; the only 
nonzero terms have oppositely charged monopoles at / and j. 
The quantity f/,j = - log G,j is in fact equal to the eff'ective 
entropic interaction between a monopole-antimonopole pair, 
and so its limiting form for large separation can be used as an 
order parameter for the confinement transition. In a confining 
phase, Uij grows without limit as the pair are separated, and 
Gij decays exponentially. Deconfinement means that Uij tends 
to a constant (with Coulomb-law corrections), so Gij has a 
nonzero limit in the Coulomb phase. At the critical point itself 
{t = 0, z - 0), the general scaling form implies 



GijiO) ~ In - rj 



(5) 



The correlation function Gij{0) can equivalently be viewed 
as the spatial distribution of two test monopoles (of opposite 
sign) inserted into a state with zero monopole density. Mea- 
suring this quantity in a numerical simulation (such as cluster 
MC 12611 ) allows the exponent y- appearing in Eq. (|2]i to be 
determined, using Eq. (|5]l. (The power law applies only for 
separations small compared to system size, so finite-size scal- 
ing is required to extract y,.) Simulations performed at z = 
therefore allow one to find <p, leading to quantitative predic- 
tions for behavior at z > 0. 

The considerations up to this point are valid for a broad 
class of transitions from the Coulomb phase. As a specific ex- 
ample, we now consider a phase transition proposed [ Il4| to 
occur in spin ice in the presence of an appropriate perturba- 
tion. This example has the advantage of being described by 
a scalar Higgs transition, allowing a quantitative comparison 
with known critical exponents. 

A nearest-neighbor (nn) model of spin ice can be written 
in the form of Eq. ([T]i, with degrees of freedom B[ - ±5 on 
the links ^ of a diamond lattice, or equivalently spins S( on 
the sites of pyrochlore constrained to local (111) axes ^. 
The nn interactions -JS( ■ Sc disfavor monopoles, giving 
z = e^^^^. One can drive a transition into a confined phase 
with, for example, an external field h{, which couples to the 
spins through a Zeeman term S = - YithfBc = 2f h{rc) ■ S(. 
The universality class of the resulting transition depends on 
the field orientation; an example of particular theoretical in- 
terest, though questionable experimental relevance, is when 
the field has a helical structure in real space (inset of Fig. |2]i. 
Explicitly, consider a field h{r) = {h cos q ■ r,h sin q ■ r, 0), 
where q - (0, 0, 2nla) is aligned with the [001] crystal direc- 
tion and a is the fee lattice constant. This perturbation favors 
a single configuration, so there is no symmetry breaking at 
T < Tc and hence no transition for z > 0. 

As argued in Ref. lll4ll . this transition is described by a 
scalar-Higgs theory, and so can be mapped to a dual transition 
in the XY universality class 12711 . Since S consists only of 
terms acting on a single link, the sum over Bt in Eq. ([T]i gives 
.S^** - 2f /(gradf 6), where / depends on h/T. The claim that 
the helical-field transition is dual to that of the XY model is 
therefore equivalent to a claim that the XY ordering transition 
is, in this case, in the same universality class as one with / 
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FIG. 2: Plot of W(T, z = 0, L) versus T/h for various system sizes 
L, showing a crossing at Tc/h =: 3.252, indicative of a continuous 
transition. (The crossing becomes increasingly sharp for larger L, 
omitted for clarity.) The inset illustrates the perturbation causing the 
transition, an applied field with a helical structure in real space, and 
the configuration of the spins for T <k Tc. 



replaced by - cos grad^ 6 1I28I1 . The absence of a transition in 
the XY model for /ixy > agrees with the observation that 
the helical-field model has no transition for z > 0. 

An alternative route to the critical theory involves mapping 
to a model of quantum bosons in 2D ill[ll2[ ]. under which the 
helical field maps to a staggered potential |14!]. The Coulomb 
phase maps to a superfluid, while the potential causes the for- 
mation of a Mott insulator at integer filling, and a transition in 
the 2 + ID XY universality class. Nonzero monopole fugacity 
causes violations of particle-number conservation; consider- 
ing the effect on the transfer matrix leads to a perturbation of 
the form zib + /?' ). A nonzero limiting value of the correla- 
tion function G,j(0) therefore signals off-diagonal long-range 
order in the condensate 1112 1. 

The XY and quantum mappings provide not only alter- 
native perspectives on the scaling analysis, but also explicit 
values for the critical exponents. Using a combination of 
MC and series expansions, Campostrini et al. ll29ll found 
V = 0.6717(l),y3 = 0.3486(1), giving cp ^ dv - p ^ 1.6665(3). 

We performed MC simulations to test the two separate pre- 
dictions resulting from our theoretical analysis: (1) the pres- 
ence of a transition in the 3D XY universality class when 
monopoles are forbidden (z = 0), and (2) scaling relations, 
such as Eq. Q, w hen z > 0. The simulations used a clus- 
ter algorithm 1261], in which loops of spins are flipped to re- 
main in the low-energy manifold, modified to allow nonzero 
monopole density. 

To determine the location of the transition and its critical 
exponents, we focus on the uniform magnetization density 
m - (16L^)"' ILe^c^ where La is the linear size of the sys- 
tem, with 16L'' spins. This has zero expectation value in both 
phases, while the quantity W{T,z,L) - L'*(|im|-), proportional 
to the "flux stiffness" in the Coulomb phase, has zero scaling 
dimension H, obeying W{T,z,L) = Q.{tL^I\zL'l'l''). 
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FIG. 3: Determination of the correlation-length critical exponent v, 
using results at z = 0. Main figure: Log-log plot of dW(T, 0, L)/dT 
atT = Tc versus system size L, fit to oc l''''. The (blue) dashed line 
shows the best fit value of the exponent v = 0.670 + 0.006, while the 
(black) solid and (purple) dash-dotted lines show fits with v fixed to 
its values for the 3D XY and Ising universality classes respectively. 
(The transition is predicted to be described by the former; the Ising 
class is displayed only for comparison.) The best fit v agrees with 
the XY class, but not Ising. Inset: Data collapse of W(T, 0, L) versus 
rL"'' using the XY exponent v = 0.6717 and Tc/h = 3.252. 



One therefore expects plots of W(T,z - 0,L) versus T for 
different values of L to cross at Tq- This is indeed confirmed 
in Fig.|2j we find Tc/h = 3.252(1), using the crossing points 
for L < 20. The slope at the crossing, dW/dT\T=Tc,z=o, is fur- 
thermore predicted to be proportional to L'^'', as confirmed 
by Fig. [3] The fitted v is consistent with the 3D XY univer- 
sality class 12 9i] . and data collapse is found (inset) using this 
value. Scaling at nonzero z is demonstrated in Fig. HI where 
W{Tc,z,L) and L^^^'''dW/dT\T=Tc ^re shown to depend on z 
and L only through zL'^^''. The inset shows monopole density 
Pni, for which a scaling form is found by taking d/d log z of 
Eq. (|2]l. In both cases, convincing data collapse is found using 
the value 0=1 .6665 appropriate to the 3D XY universality 
class. Larger system sizes would be required to place stringent 
bounds on the crossover exponent; we estimate (p = 1.65(15), 
with confidence interval based on the quality of data collapse. 

With scaling at nonzero monopole fugacity established in 
the helical-field transition, we now consider prospects for its 
demonstration in other models and in experiment. MC re- 
sults for spin ice in a (100) field have been reported previously 
I lin . including at nonzero monopole fugacity. The transition 
is at its upper critical dimension, so logarithmic corrections 
to expressions such as Eq. ^ are expected; an extension of 
the current theory and comparison with numerics will be pre- 
sented elsewhere. Our predictions also apply to the classi- 
cal dimer model, where monopoles correspond to empty or 
doubly-occupied sites, and can be tested in numerical simula- 
tions ISllaLZll extended to allow such defects. 

The leading candidates for realizing this class of transition 
are the spin ices ||9lll0|], such as Dy2Ti207 (DTO), in which 
the dipolar correlations 113 Oil and monopole excitations [Ift] of 



FIG. 4: Scaling at nonzero monopole fugacity z- Main figure: Plot 
of W(Tc,z, L) (empty symbols, left scale) and dW{T, z, L)/dT atT = 
Tc (filled symbols, right scale) versus z^L"^- The data collapse onto a 
single curve in each case, with v and (p = 1.6665 taking their values 
for the 3D XY universality class. Inset: Monopole density pni divided 
by |«p"", using the 3D XY exponent a = -0.015 (H]. Plotted against 
z|/r**, data for various T and z (and fixed L = 16) collapse separately 
for each sign of t. 



the Coulomb phase are well established. Dynamical freez- 
ing at a temperature Tf ^ 0.6 K sets a lower limit of roughly 
z ^ 10"^ (using A ^ 4K ^^ for DTO), suggesting the pa- 
rameter regime of Fig. |4] may be accessible. While these 
simulations apply to the experimentally challenging case of 
a helical field, it seems reasonable to hope that, since one can 
achieve z low enough to observe the Coulomb phase, clear 
signatures of a confinement transition may also be visible in, 
for example, spin ice in a (100) field l32ll . If so, the "smok- 
ing gun" for this type of transition would be the distinctive 
scaling with z = e"^^^ of, for example, the heat capacity or 
the monopole density p,^ (which determines the broadening 
of "pinch points" in neutron scattering). 

The analysis presented here is directly applicable to quan- 
tum systems, such as the "quantum spin ice" Yb2Ti207 
1 19l|33|], in the case of a nonzero- temperature continuous tran- 
sition from the Coulomb phase. A similar approach may be 
applicable to quantum phase transitions in 2D, where the U(l) 
spin liquid is also unstable to monopoles JITII . 

In summary, this work has provided a unified perspective 
on a family of unconventional transitions out of the Coulomb 
phase, incorporating the effects of monopole defects on the 
critical behavior The main contributions are: (1) a general 
theory of universal scaling at nonzero monopole fugacity in 
terms of the properties of the unconventional fixed point, (2) 
a detailed understanding of these phenomena for a particular 
transition, that of spin ice in the presence of a helical magnetic 
field, using mappings to the XY model and quantum bosons, 
and (3) confirmation from numerical simulations that the pre- 
dicted helical-field transition exists, that it obeys scaling rela- 
tions both without and with monopoles, and that the exponents 
are consistent with those of the 3D XY universality class. 
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